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Abstract. Doppler CORALIE measurements of the solar-type stars HD 141937, HD 162020, HD 168443 and 
HD 202206 show Keplerian radial- velocity variations revealing the presence of 4 new companions with mini- 
mum masses close to the planet/brown-dwarf transition, namely with rri2 sin i = 9.7, 14.4, 16.9, and 17.5 Mj up , 
respectively. The orbits present fairly large eccentricities (0.22 < e < 0.43). Except for HD 162020, the parent 
stars are metal rich compared to the Sun, as are most of the detected extra-solar planet hosts. Considerations 
of tidal dissipation in the short-period HD 162020 system points towards a brown-dwarf nature for the low-mass 
companion. HD 168443 is a multiple system with two low-mass companions being either brown dwarfs or formed 
simultaneously in the protoplanetary disks as superplanets. For HD 202206, the radial velocities show an addi- 
tional drift revealing a further outer companion, the nature of which is still unknown. Finally, the stellar-host and 
orbital properties of massive planets are examined in comparison to lighter exoplanets. Observed trends include 
the need of metal-rich stars to form massive exoplanets and the lack of short periods for massive planets. If 
confirmed with improved statistics, these features may provide constraints for the migration scenario. 
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1. Introduction 

Since the discovery of the extra-solar planet orbiting 
51 Peg (Mayor k Queloz 1995), high-precision radial- 



velocity measurements proved to be very efficient for de- 
tecting very low-mass companions to solar-type stars. In 
about 6 years, close to 80 planetary candidates with min- 
imum mass 777,2 sin i < 10Mj up have been announced in- 
cluding 7 planetary systems and a few sub-Saturnian plan- 
ets (see e.g. [Udry k Mayor 200 1|; [Fischer et al. 2002j, for 



recent reference updates). 

Interestingly, brown-dwarf candidates, easier to detect 
with high-precision Doppler surveys, seem to be more 
sparse than exoplanets (Mayor et al. 1997), especially in 
the 10-30M Jup interval ( palbwachs et al. 2000b| ), the so- 
called brown-dwarf desert. Objects in this domain are very 
important to understand the brown-dwarf/planet transi- 
tion. The paradigm behind the distinction between plan- 
ets and brown dwarfs may rely on different consider- 
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ations: mass, physics of the interior, formation mecha- 
nism, etc. From the "formation" point of view, the brown- 
dwarf companions belong to the low-mass end of the sec- 
ondaries formed in binary stars whereas planets form 
in the protostellar disk. Such distinct origins of plane- 
tary and multiple-star systems is clearly emphasized by 
the two peaks in the observed distribution of minimum 
masses of secondaries to solar-type stars as shown in Fig. [I] 
(top), providing an updatec^] version of the diagram. They 
strongly suggest different formation and evolution histo- 
ries for the two populations. Below 10Mj up the planetary 
distribution increases with decreasing mass and is thus not 
the tail of the stellar binary distribution. 

For the objects detected by the radial-velocity tech- 
nique, only minimum masses are determined because the 
inclination angles of the orbital planes relative to the line 
of sight cannot be derived from the spectroscopic data 
only. The determination of the true masses for most of 
these objects is expected soon with high-angular resolu- 



All candidates known on the 15th of November 2001. The 
masses are from the discovery papers or from CORALIE- 
ELODIE orbital solutions. A summary table is provided at 
obswww. unige. ch/~naef/who_discovered_that_planet. html 



S. Udry et al.: The CORALIE survey for southern extra-solar planets VIII. 



0.6 



-OA 

0) 

a 
cr 

0! 




2000a). The CORALIE survey has been very successful, 



with the detection of a significant fraction of the known ex- 
oplanet candidates (Udry & Mayor 2001). Information on 
the method, technical and instrumental details are given 
in Baranne et al. (1996| ). Recent improvements in the re- 



duction software have allowed us to bring the long-term in- 
strumental precision of individual measurements from the 
previously obtained ~ 7 ms _1 down to ^2 ms _1 (Queloz 



et al. 2001b). Asteroseismology measurements of aCenA 
even show a short-term precision below the 1 ms _1 limit 
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over 1 night (Bouchy & Carrier 2001). The size of the tele- 
scope is now the main limitation of the precision actually 
achieved for most of our sample stars. 

The cross-correlation technique used has proven to be 
very simple, robust and efficient for radial-velocity mea- 
surements. However, it is still not optimum in terms of 



Fig. 1. Top: Observed distribution of minimum masses of 
secondaries to solar-type stars (log scale). Stellar binaries 
with G and K primaries are from Duqucnnoy & Mayor 
(19910 and [Halbwachs ct al. (2000a| ). Bottom: Updated 
statistical distribution of true giant-planet masses (dashed 
line; derived as in lorissen et al. 2001) superimposed on 
the ?7i2 sin i planetary distribution 



tion astrometric facilities that will become available in 
the upcoming years (e.g. PRIMA on the VLTI, SIM). It 
is however already possible to apply a statistical decon- 
volution to the growing sample of exoplanet candidate s 
fljorisscn, Mayor, fc Udry 200 1| ; |Zuckcr fc Mazeh 200 ip . 
The updated planetary true mass distribution, derived as 
in Jorissen et al. (2001 ), is presented in Fig. [l] (bottom). 
Probably because of the strong observational bias favour- 
ing the more massive planets, the distribution looks "bi- 
modal" . A careful treatment of the bias, however, has to be 
done before being able to convincingly interpret the shape 
of the distribution. Nevertheless, the bias at the high-mass 
end of the planetary distribution is vanishing and the ob- 
servational maximum mass of exoplanets is fairly well de- 
termined, around 10-11 Mj up where the curve drops to 
almost zero. If this value is correct then the question of 
the true masses and nature of the candidates with slightly 
higher values of rr^sini is becoming very interesting. 

Since the summer of 1998, a large high-precision 
radial-velocity programme has been carried out with the 
CORALIE echelle spectrograph on the 1.2-m Euler Swiss 



Doppler information extraction from the spectra (Bouchy 
ct al. 200 1| ; phelli 2000[ ). Recently, |Pepe et al. (2002| ) im- 
proved the procedure by introducing a "correct" weighting 
of the spectral lines involved in the cross correlation. At 
the same time they also reduced the astroclimatic-induced 
noise by restricting more severely the zone of the spectra 
potentially affected by telluric lines. The overall gain cor- 
responds to a virtual increase in the signal-to-noise by a 
factor of ~ 1.25 (i.e. a virtual decrease of the photon- noise 
error by the same factor). 

This new procedure for radial-velocity estimation is 
used for the objects described in this paper. We present 
4 new very low-mass companions to solar-type stars, de- 
tected with CORALIE, and with minimum masses in the 
planet-brown dwarf transition domain. With four other 
objects - HP 114762 ( [Latham et al. 1989|), HP 11 0833^ 
([Mayor et al. 1997|), HP 39091 ([Hugh ct al. 2002|) and 



HP 136118 ( Fischer et al. 2002| ) - these candidates are 
the only companions to solar-type stars known to date 
with minimum masses between ^10 and 20Mj up . They 
are thus of prime importance for the description of the 
transition zone between brown dwarfs and planets. The 
first sections of the paper are dedicated to the descrip- 
tion of the stellar properties of the hosts of the the new 
candidates, then to their orbital characteristics and to a 
discussion of the possible nature of these objects. Finally, 
the orbital properties and the characteristics of stars with 
"massive" planets are examined in comparison with sys- 
tems harbouring "lighter" planets. 

The radial- velocity data for the new candidates will be 
made available in electronic form at the Centre de Donnees 
Stellaires (CDS) in Strasbourg. 

2. Stellar characteristics of the candidate hosts 

The stars hosting the 4 very low-mass candidates pre- 
sented here were observed by the HIPPARCOS astrometric 
satellite. Most of the quoted photometric and astrometric 
parameters are thus taken from the mission output cata- 



HP 110833 was shown to be a stellar binary by Halbwach 



telescope at LaSilla ( |Queloz et al. 2000b[ ; |Udry et al] |et al. (20()0h] ) using the HIPPARCOS astrometric data 
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Table 1. Observed and inferred stellar parameters for HP 14 1937, HD 162020, HD 168443 and HD 202206. Photometric 
and astrometric parameters are from HIPPARCOS (ESA 1997). The atmospheric parameters T c g, log g, [Fe/H] are from 
Santos et al. (200"la| ,p|) and Gonzalez et al. (2001). Spectral types are from HIPPARCOS or derived from spectroscopy. 
The bolometric correction is computed from Flower (1996) a using the spectroscopic T c g determinations. The projected 
rotational velocities come from a calibration of the CORALIE cross-correlation functions. Activity indicators (Scor, 
\ogR' HK ) are estimated from N act high S/N CORALIE spectra, following Santos et al. (2000| ). The given ages are 
derived fro m the \ogR' HK activ ity indicator (Donahue 1993) or/and from the Geneva evolutionary models (Schallci 
et al. 1992; Schaerer et al. 1993) which also provide mass estimates. Finally, rotational periods are obtained from the 
logi?^^ as well, following Noyes et al. (1984) 



Parameter 




HD 141937 


HD 162020 


HD 168443 


HD 202206 


Spectral Type 




G2/G3V 


K3V 


G8IV" 


G6V 


V 




7.25 


9.10 


6.92 


8.08 


B-V 




0.628 


0.964 


0.724 


0.714 


TV 


[mas] 


29.89 ± 1.08 


31.99 ± 1.48 


26.40 ± 0.85 


21.58 ± 1.14 


M v 




4.63 


6.63 


4.03 


4.75 


BC 




-0.055 


-0.388 


-0.125 


-0.082 


L 


[Le] 


1.17 


0.25 


2.17 


1.07 


[Fe/H] 




0.11 


0.01 ±0.11 


0.10 ± 0.03 


0.37 ±0.07 


M 


[M Q ] 


1.1 


0.75 


1.01 6 


1.15 


T c s 


[K] 


5925 


4830 ± 80 


5555 ± 40 


5765 ± 40 


logs 


[cgs] 


4.62 


4.76 ±0.25 


4.10 ± 0.12 


4.75 ±0.20 


v sini 


[kms" 1 ] 


2.1 


1.9 


1.7 


2.5 


iVact 




19 




1 




Scon 




0.23 ±0.04 




0.19 




Smw 




0.24 




0.21 




log(-RW) 




-4.65 




-4.8 (-5.08 6 ) 




Piot(RHK) 


[days] 


13.25 




26.8 (37 6 ) 




age (R' HK /Model) 


[Gyr] 


1.6/2 ±6 


-/6 ± 18 


7.8 b /> 10 


-/5.6±1.2 



"Quoted values i n the paper include errors. The correct values have been obtained directly from the author 
'Value quoted in |Marcy et al. (1999|) 



logue (ESA 1997). High-precision spectroscopic studies of G2V spectral type and color index [B—V = 0.628). SantO: 



these stars have also been performed by several authors et al. (2001b) performed a high-resolution spectroscopic 



in the context of examining the metallicity distribution 
of stars hosting planets in comparison to "single" stars 
of the solar neighbourhood ( Gonzalez et al. 2001 ; |Santos 



et al. a 001a| . b|) . Observed and inferred stellar parameters 



from these different sources are summarized in Table [jj. In 
the table, the given masses of primary stars are estimated 
from evolutionary tracks of the Geneva models with ap- 
propriate spectrosco pic parameters (TeS , L, metallicity; 
Schallcr ct al. 1992 ; |Schaerer et al. 1993| ). Age estimates 
are also provided by those models but usually with very 
large uncertainties for our type of stars. The projected 
rotational velocity, wsini, comes from the calibration of 
the CORALIE cross-correlation functions^ (CCF) derived 
in the same way as the calibration of the ELODIE CCF 
dQucloz ct al. 1998| ). 



2.1. HD 141937 (HIP 77740) 

From the HIPPARCOS parallax (29.89 ± 1.08 mas) and vi- 
sual magnitude (V — 7.25), we derive for HD 141937 an ab- 
solute magnitude My — 4.63 in agreement with its given 



The calibration does not account for metallicity effects. For 
metal-rich stars the rotational broadening is therefore slightly 
overestimated 



abundance study for this star and derived precise values 
for its effective temperature (T e ff = 5925 K), metallicity 
([Fe/H] = 0.11) and gravity (log<? = 4.62), using a stan- 
dard local thermodynamical equilibrium (LTE) analysis. 
Using a calib rated bolometric correction BC = —0.055 
( Flower 1996 ) combined with the spectroscopic T e g deter- 
mination, the star luminosity is found to be L = 1.17Lq. 
A mass M — 1.1 M Q and a badly constrained age of 
2 ± 6 Gyr are then estimated from the Geneva evolution- 



ary models (Bchaller et al. 1992). 

The dispersion of the HIPPARCOS photometric data 
(effp = 0.007 mag) shows no evidence of variation of the 
star luminosity at the instrument precision. The same con- 
clusion holds from the Geneva photometry observations. 

2.2. HD 162020 (HIP 87330) 

In the HIPPARCOS catalogue HD 162020 is a K2 dwarf 
with V = 6.35 and B - V = 0.964. The catalogue also 
lists a precise astrometric parallax it = 31.99 ± 1.48 mas 
corresponding to a distance of 31.26 pc from the Sun. 
The derived absolute magnitude My — 6.63 is typical 
for a K3 dwarf. From high-resolution CORALIE spectra, 
pantos et al. (2001a ) derived the following spectroscopic 
parameters: T cff = 4830 K, [Fe/H] = 0.01 and log 5 = 4.76 
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(Table [l]) . With a calibrated bolometric correction BC = 
—0.388 ( [Flower 1996| ) and the derived effective tempera- 
ture the star luminosity is found to be L — 0.25 L , also 
suggesting a K3V spectral type. Solar-metallicity models 
( (Bchaller et al. 1992j ) give then a mass of 0.75 M . They 
also point towards an old stellar age but with a very large 
uncertainty, much larger than the age itself (6±18Gyr). 
On the other hand, we will see further that the activity 
level of the star rather suggests a younger age. 

The dispersion of the HIPPARCOS photometric data of 
HD 162020 {a Hp = 0.018 mag) is found to be slightly higher 
than the value expected for a 9th-magnitude star. We will 
see below that the star shows spectral indications of ac- 
tivity that can explain this feature. 



2.3. HD 168443 (HIP 89844, G J 4052) 



The stellar characteristics of HD 168443 have been dis- 



cussed by Marcy et al. (199£ ) in the paper announcing 
the detection of the first planet orbiting the star. From 
the HIPPARCOS photometric data, they have estimated 
the star to reside ~ 1.5 mag above the main sequence 
in the HR diagram (My — 4.03), being slightly evolved. 
They estimate the star to have a G8IV spectral type and 
a mass of 1.01 M Q . A recent high-precision spectroscopic 
study by ponzalez et al. (2001 ) completely corroborates 
that result by deriving reliable values of the effective tem- 
perature (T off = 5555 K), metallicity ([Fe /H]=0.1) and 
gravity (log 17 = 4.1) of the star. Using Flower (1996) 
calibration for the bolometric correction (BC = —0.125), 
one obtains a luminosity L = 2.17L . The mass derived 



from the Geneva evolutionary tracks (Schallcr et al. 1992) 
agrees with the estimate of Marcy et al. (1999] ) and the age 
is found to be larger than lOGyr, suitable for a slightly 
evolved late G star. These main stellar properties are re- 
called in Table @. 

The star is found to be photometrically stable in both 
the HIPPARCOS data (a Hp = 0.007 mag) and the Geneva 
photometry (ay = 0.005 mag). 

2. 4. HD 202206 (HIP 1 04903 ) 



In the HIPPARCOS catalogue, HD 202206 is given as a 
G6 dwarf of visual magnitude V = 8.08 and color index 
B — V = 0.714. The measured parallax (21.58 ± 1.14mas) 
leads to an absolute magnitude My — 4.75, ~ 0.4 mag 
brighter than the expected value for a typical G6 dwarf 



of solar metallicity. From CORALIE spectra, Santos et al 



(2001a) derived an effective temperature T e g = 5765 K, 
a gravity logg = 4.75 and a very high metal content 
[Fe/H]=0.37 (Table 0). The very high metallicity of 
HD 202206 probably accounts for its overluminosity as T e s 
is also larger than the value expected for a G6 dwarf. From 
BC = -0.082, we derive L = 1.07 L and the Geneva 
models (Bchaerer et al. 1993) yields a mass of M = 1.15 M 



As for HD 162020, the dispersion of the HIPPARCOS 
photometric data of HD 202206 (a Hp = 0.013 mag) is a bit 
high for the star magnitude but again some indication of 
stellar activity is seen in the spectra (see below). 

2.5. Chromospheric activity 

The amplitude of the radial- velocity jitter associated with 
intrinsic stellar activity may reach a few tens of ms _1 , 
especially for high-rotation stars with large spectral-line 
asymmetry due to spots ( fjaar fc Donahue 1997 ). The 



coherent-spot survival on the stellar surface over several 
rotational periods may even mimic the radial- velocity vari- 



ation induced by a planetary companion (Qucloz et al 
|2001a| ). 



The stellar activity can be associated with the pres- 
ence of chromospheric emission in the centre of the Ca n H 
and K absorption lines. When reported to the photo- 
spheric flux, the intensity of this emission provides good 
quantitative estimators (S'corQ or log^HK) °f tne ac ~ 
tivity level of solar-type stars and thus to the expected 
level of induced spurious noise on the radial- velocity mea- 
surements (Saar et al. 1998; Santos et al. 2000). For the 
brightest stars in our planet-search sample, such indicators 
are directly measured on the CORALIE spectra. However, 
the low brightness of some of the targets does not al- 
low us to derive good indicators^ (a s e.g. for HD 162020 
and HD 202206; pantos et al. 2000[ ). In such cases, the 
available spectra are added to obtain a high signal in the 
A 3968.5 A CailH absorption line region. The resulting 
spectrum may contain traces of the Thorium- Argon spec- 



trum used as a radial- velocity reference (Baranne et al 



1996) and is thus not optimal for precise spectroscopic 



studies in this region but it at least allows for a visual 
check of the chromospheric emission in the center of the 
line. Figure^ shows the corresponding spectral domain for 
the stars presented in this paper. 

Moderate values of the \ogR^ K chromospheric activ- 
ity indicator have been derived for HD 141937 (—4.65) and 
HD 168443(-4.80) (Table g see |Santos et al. 200C| , for 
details about the technique). Moreover, the correspond- 
ing A 3968.5 A CailH absorption lines do not show clear 
chromospheric emission features (Fig. ||). Radial- velocity 
jitters are thus not expected to be large for these 2 stars al- 
though the effect might be slightly increased for the early 
G dwarf HD 141937 showing a non-zero projected rota- 
tional velocity (t>sini = 2.1 kms -1 ). 



4 The 5*cor activity index has been calibrated t o the Mount- 
Wilson system index Smw ( Vaughan et al. 1978 ) to compute 



and an age of about 5.6 ±1.2 Gyr. 



the log -Rhk indicator ( Noyes et al. 1984^ 

The correction for the diffuse light in the spectrograph is 
not known precisely enough and for low S/N spectra it rep- 
resents a significant fraction of the light in the line. It then 
induces systematic errors in the estimate of the activity indi- 
cators 

6 Our esti mate only rests on 1 high-S/N spectrum. Marcy 
^t al. (1999|) quote a more reliable value of —5.08 suggesting 



that HD 168443 is an even more quiet star 
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3965 3970 
\(A) [Ca II H] 

Fig. 2. A 3968.5 A CaiiH absorption line region of the 
summed CORALIE spectra for the 4 stars considered 
in this paper. Clear emission features are observed for 
HD 162020 and HD 202206. For the latter a small part of 
the HD 168443 spectrum in the center of the line has also 
been reported in the diagram (dotted line) to emphasize 
its chromospheric emission. For clarity, spurious emission 
contamination features of the Thorium- Argon lamp have 
been removed from the spectrum of HD 162020. No trace of 
chromospheric emission is visually observed for HD 141937 
and HD 168443, in agreement with their moderate mea- 
sured values of log i?^ K (Table |l|) 



Although noticeable, the activity level of HD 202206 is 
not very important (Fig. ^) and should not cause any trou- 
ble beyond adding some low-level high-frequency spurious 
noise in the radial- velocity measurements. 

Finally, from the log- R^ K value (when a vailable), fol- 
lowing the calibration by Noyes et al. (1984 ) we can esti- 
mate the rotational period o f the star as wel l as the stellar 
ag e using the calibra tion in Donahue (1993 ) (also quoted 
in Henry et al. 1996] ). The inferred statistical values are 
given m Table |l|. The so-derived ages are compatible with 
values provided by evolutionary tracks, taking into ac- 
count the large uncertainties in the age determinations. 

3. HD 141937 orbital solution 

Over 882 days, we obtained 81 CORALIE observations of 
HD 141937, with a photon-noise uncertainty distribution 
peaking around 6 ms _1 . A fairly large long-term radial- 
velocity variation was noticed early on, but we had to wait 
for more than one orbital period to derive good orbital pa- 
rameters because the star was unfortunately behind the 
Sun at the time of the maximum and minimum of the 
radial-velocity curve. The best Keplerian fit to the data 
yields a precise period P = 653.2 days, an eccentricity 
e = 0.41 and a semi-amplitude if = 234.5 ms" 1 (Table ||). 
The radial-velocity curve is displayed in Fig. ^ with the 
residuals around the solution. The mass M — 1.1 M© de- 
rived for the parent star leads to a planet minimum mass 
7712 sini = 9.7Mj up . 

The weighted r.m.s. around the best Keplerian solution 
(er(0 — C) = 8.7 ms -1 ) is large compared to the typical 
measurement uncertainties (x? e d = 2.69). As no longer- 
term variation clearly arises from the residuals of the fit 
(Fig. ||), the reason for the slight extra noise (~ 6ms" 1 ) 
probably lies in the star activity (log R' HK — —4.65) cou- 
pled with the observed non-zero stellar rotatio n (v sin i = 
2.1 km s~ 1 ) , in agreement wit h values quoted by [5aar et al 
Kl998| ) or [Santos et al. (200C| ). 



Although too faint to provide a reliable estimate of 
\ogR' HK , HD 162020 clearly shows a strong emission fea- 
ture in the core of the CaiiH line. The activity-related 
radial- velocity jitter of a slowly rotating K2 dwarf is how- 
ever not expected to be large ( Santos et al. 200C ). It is 
certainly not responsible for the large radial- velocity varia- 
tion observed for HD 162020 (3.3 kms" 1 peak to peak). As 
shown by Henry (199S), activity could on the other hand 
be invoked to explain the dispersion of the HIPPARCOS 
photometric data. The activity level points towards a 
young age for the star. However, we will see in Sect. 4.1.1 
that the system has probably been synchronized over the 
lifetime of the star, thus increasing the stellar rotation. 
As activity is very sensitive to rotation, the activity level 
could also have been boosted by the synchronization mak- 
ing the star look younger than its real age. 



4. HD 162020 b: A "Hot Brown Dwarf" 

Between the 24th of June 1999 and the 14th of October 
2001, 46 CORALIE radial velocities of HD 162020 were 
gathered. The low brightness of the star limits the photon 
noise of our measurements to about 8 ms -1 in a typi- 
cal integration time of 15 minutes, under normal weather 
conditions. As the observed radial-velocity variation was 
large, several observations were made under worse me- 
teorological conditions and the distribution of measure- 
ment errors shows a right-end tail up to 40 ms" 1 . The 
short orbital period and large observed semi-amplitude of 
the radial-velocity variation allowed us, however, to very 
rapidly^ determine orbital elements and characteristics of 
the companion. The imprecise subset of measurements 
limit the quality of the solution (r.m.s. ~ 13.6 ms" 1 , 



7 The discovery was announced on the 4th of May 2000 by 
an ESO press release (|ESO 2000[) 
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Fig. 3. Top. CORALIE radial-velocity measurements of 
HD 141937 superimposed on the best Keplerian model. 
Error bars (photon noise) are very small in compari- 
son with the amplitude of the radial-velocity variation. 
Bottom. Residuals around the solution 



Xrcd = 2.46), and we have derived a new solution with 
only the 30 CORALIE observations with photon-noise er- 
rors below 10 ms _1 . This solution yields a short period 
P of 8.428198 ± 0.000056 days with a non-zero eccentric- 
ity e = 0.277 ± 0.002. Taking 0.75 M as the primary mass 
(see above), the derived orbital parameters lead to a min- 
imum mass n^sini = 14.4 Mj up for the companion. The 
very short orbital period derived for HD 162020 locates the 
companion only 0.074 AU from the primary star. At such 
a small distance from its parent star and following e.g. 
Guillot ct al. (1996), the companion equilibrium temper- 

The 



ature at the surface is estimated to be around 650 K. 
companion is thus a "hot" superplanet or brown dwarf on 
a non-circular orbit. The complete set of orbital elements 
with their uncertainties is given in Table ^ as well as some 
interesting inferred quantities. 

The weighted r.m.s. to the Keplerian fit is 8.1ms _1 
and the reduced x 2 of the solution is 1.37. The improve- 
ment from the preliminary solution including all mea- 
sureme nts clearly shows the significant degrading effect of 



Fig. 4. Top. Phased CORALIE radial-velocity measure- 
ments and Keplerian orbital solution for HD 162020. The 
solution is calculated with only the 30 higher signal-to- 
noise measurements but all the 46 radial velocities with 
photon-noise errors (error bars) are reported on the dia- 
gram. Bottom. Residuals around the solution displayed as 
a function of time 

4.1. Tidal dissipation in the HD 162020 system 

In the same way as close binaries, giant gaseous plan- 
ets closely orbiting their stars are subject to spin-orbit 
synchronization and/or orbital circularization associated 
with tidal dissipation in the star or the planet. The tilted 
mass distribution induced in the convective envelope of 
the object by the gravitational attraction of the compan- 
ion is phase shifted because of dissipation. It then exerts 
a torque on the companion, leading to an exchange of an- 
gular momentum between its spin and the orbital motion, 
tending to synchronize and circularize the orbit. When 
the tidally distorted star has a convective envelope, the 
tidal dissipation may be represented by the viscosity of 
convective eddies ( ^ahn 1989| ). 

For binaries with dwarf-star primaries, orbits with pe- 
riods smaller than ~ 10 days are circularized (see e.g. 
Mayor et al. 2001, for a recent review). Much attention 



the lo wer-quality velocities on the derived solution. The 

residuals around the solution show no significant evidence 
of a possible long-period additional companion (Fig. |J). 
Although the star shows clear indication of chromospheric 
activity, the level of activity-induced extra noise is very 
small, in agreement with values predicted for slowly ro- 
tating K dwarfs (3aar et al. 199S; Santos et al. 2000). 



has also been paid to this eff e ct for hot Jupitc rs (Rasic 
ct al. 1996| ; |Marcy et al. 1997| ; [Ford et al. 1999| ). We will 



focus here on the case of HD 162020 in which the "planet" 
is slightly more distant than the one orbiting 51 Peg but 
also much heavier. 

4.1.1. Tidal dissipation in the star and synchronization 



Following the treatment of Zahn (1989, 1992) for objects 
with a convective envelope, in the case of not too high 
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eccentricities, the time scale for synchronization by tidal 
dissipation in the star through convective viscosity is given 

by 



l 



A 2 2 M+Rl (R* 
tf I* \ a 



(1) 



where q = m p \/M i < is the planet to star mass ratio, i?* 
and 1+ are the stellar radius and moment of inertia, a is 
the orbital separation, tf is the characteristic time of the 
physical process responsible for the dissipation (turnover 
timescale of the eddies) and A2 is close to the apsidal con- 
stant measuring the response to the external torque im- 
posed by the companion when the star is fully convective 
and smaller when the star has a finite convective zone. 
Zahn (199^ ) has tabulated fc 2 = J*/(M*fl 2 ), tf and A 2 



for the zero-age main sequence. Interpolating in his table, 
with M* = O.75M , we obtain for HD 162020: fc 2 = 0.132, 
£/ = 0.509yr and A 2 = 0.0094. The ty pical time t f is usu- 
ally defined as t f = (M+Rl/L*) 1 / 3 ( fcahn 1994|) . In the 
case of non-fully convective dwarfs, an estimate of 



tf — (M env R env (R+ i? entJ )/3Z/^) ^ , 



(2) 



where M env is the mass of the convective envelop e and 
Renv t he ra dius at its base, is more appropriate ( Rasio 
et al. t996| ). For HD 162020 {M env /M+ ~ 0.17, |Murray 
et al. joOlj, Renv/R* - 0.7), we calculate t' f ~ 0.13yr. 



With q = 0.01834/ sinz, we derive 
Tsync (hd 162020) ~ 1.3 • 10 10 sin 2 i yr. 



(3) 



Even taking into account our limited knowledge of the 
viscous dissipative process and the large uncertainty of the 
stellar age determination, it appears that HD 162020 can 
be synchronized by the close-orbiting low-mass compan- 
ion. In such a case, when the two components are close, 
t he synchr onized state may be unstable to orbital decay 
( Hut 1980 ). Instability occurs when the ratio of spin to or- 
bital angular momentum J sp m/</orb > 1/3. For HD 162020, 
•/spin/^orb — 0.007. Hence, the star was brought into syn- 
chronicity without reaching tidal instability. 

In the case of synchronization, knowing the stellar ro- 
tational period and assuming that the orbital and rotation 
axes coincide^], the equatorial velocity V cq may be derived 
from the star radius and then the orbital plane inclination 
is obtained from the measured projected rotational veloc- 
ity. Using a typical value i?+ ~ 0.75 Rq for a K3 dwarf, 
we get V cq ~ 4.5 kms" 1 and an indicative sini ~ 0.42 
(usini = 1.9 kms -1 ). This leads to m 2 — 34M,j up i.e. a 
probable brown dwarf for the companion of HD 162020. Of 
course, the uncertainty on the usini value is difficult to 
estimate and probably it does not exclude the companion 
from having a low stellar mass. 



As observed e.g. for HD 209458 (Dueloz et al. 2000a) 



However, it should be pointed out that, in the case of a brown- 
dwarf companion, this may not be the case 



4.1.2. Circularization due to tidal dissipation in the star 

The observed significant orbital eccentricity of HD 162020 
(e = 0.277 ± 0.002) shows that the circularization induced 
by tidal dissipation in the stellar convective envelope had 
no time to proceed over the age of the star. Assuming 
again standard tidal dissi pation theory (e.g. Zahn 1989| ; 



Verbunt fc Phinncy 1995), an explicit expression for the 



circularization time is given by 
dine / M env 1 



1 



dt 



t' f M* q(l + q) 



a 



leading to 

T c i r (HD162020) 



4.5 -10 9 ' 



/ 



yr. 



(4) 



(5) 



The parameters are the same as above and / is obtained 
by integrating viscous dissipation of tidal energy through- 
out the convective zone. / ~ 1 as long as t' f << P or b 



(e.g. ^ahn & Bouchct 1989 


; Verbunt & Phinney 1995 


Rasio et al. 1996; 


Ford et al. 1999 


) corresponding to the 



case where the main contribution to the viscosity comes 
from the largest convective cells. As is clearly explained 
in Rasio et al. (1996] ), for P or b < t' f , the largest eddies can 



no longer contribute to the viscosity because the velocity 
field they are damping will have changed direction before 
they can transfer momentum. Only eddies with turnover 
times smaller than the pumping period (P or b/2) will con- 
tribute and the eddy viscosity is then reduced by a fac- 
tor {2t'f/P or b) a . The value of a is d ebated but generally 
thoug ht to be 1 (|Zahn 1992|) or 2 ( |Goldrcich k Keeley 
1977). So, in general we have 



/ = / 'min[l,(^r 



(6) 



with f oil. Choosing a — 2 (the extreme case), we esti- 
mate then for HD 162020 / ~ 0.015 and the circularization 
time becomes 



(HD 162020) 



10 11 sini yr. 



(7) 



Even taking into account the above rough estimate of 
sini, the tidal dissipation in the convective envelope of 
HD 162020 is thus not supposed to have circularized the 
orbit. On the other hand, the derived t c - it value does not 
allow for a very small value of the sinz. 

4.1.3. Tidal dissipation in the planet 

Can tides in the low-mass companion have been more ef- 
ficient than stellar tides in circularizing the orbit, as it is 
the case for 51 Peg ( Rasio et al. 1996 )? In this case, the 
typical circularization time is given by 



4 

63 ( 



GM+ 



1/2 



/// 



(8) 



(Rasio et al. 1996, and quoted references), where Q is 
proportional to the tidal pumping period (P or b/2) and is 
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Table 2. CORALIE best Keplerian orbital solutions derived for HD 141937, HD 162020 and HD 202206 as well as 
inferred planetary parameters. Note that the parameter uncertainties are directly taken from the diagonal elements of 
the covariance matrix as if the parameters were uncorrelated. These uncertainties are thus probably underestimated 



Parameter 


x x j_/ initio i 


nri i R2n2n 

111/ J U Zi\J Zi\J 


111./ ^u^^uu 


P 


[days] 


653 22 + 1 21 


8 4281 98 + 000056 


256 003 + 062 


T 


[JD- 2400000] 


51847 38 ± 1 97 


51990 677 ± 005 


51919.02 ±0.16 


e 




41 +0 01 


V ' . — 1 I _1_ U.UUi 


429 + 002 


V 


[kms" 1 ] 


—2.915 ± 0.002 


—27.328 ± 0.002 


14.681 ±0.002 


w 


[deg] 


1 87 72 + 80 


28 40 + 23 


160 32 +0 31 


K 


[ms^ 1 ] 


234.5 ± 6.4 


1813 ± 4 


564 8 ± 1 3 


Linear drift [ms x yr 1 ] 






42.9 ±1.3 


N meas 




81 


30 


95 


cr(0 - C) 

Xrcd 


[ms" 1 ] 


8.7 


8.1 


9.5 




2.69 


1.37 


2.38 


ai sin i 


[AU] 


0.01282 


0.00135 


0.012 


f(m) 


[10- 6 M o ] 


0.658 


4.620 


3.526 


m2 sini 


[M Jup ] 


9.7 


14.4 


17.5 


a 


[AU] 


1.52 


0.074 


0.83 


7eq 


[K] 




650 





about 10 5 for Jupiter. Comparing 51 Peg and HD 162020 
we can write 



r 51 Peg 
362020 



a 



51 Pee 



p51 Peg 

orb 

— pl62020 I „162020 

r orb \ a 

51 Peg / n 

x m pl t I R 

x ^,162020 



13/2 



^162020 



3/2 




(9) 



From the value quoted in Rasio ct al. (1996 ) we obtain 



N „ „ ia sini(162020) 

T e (HD 162020) ~ 2 • lO 1 ^ . , ( JT. 

smi(51 Peg) 



(10) 



The tides in the low-mass companion were thus not able 
to circularize the orbit over the lifetime of the star which 
is expected from the observed non-zero value of the orbital 
eccentricity. 

In conclusion, it appears that the companion of 
HD 162020 is probably a brown dwarf (although a low- 
mass star cannot be ruled out). The observed eccentricity 
is so not surprising if the system was formed as for usual 
binary stars. This system will be a perfect target for fu- 
ture high-angular resolution astrometric facilities. It po- 
tentially will provide the true mass of a low-mass brown 
dwarf. 

4.2. Photometric transit search 

The short orbital period makes HD 162020 a good target 
for a photometric transit search. The star was intensively 
followed with the Danish SAT at LaSilla in collabora- 
tion with our colleagues from Copenhagen, unfortunately 
without success. This is not a surprise in view of the con- 
siderations developed in the previous sections about the 
probable orbital inclination of the system. The photomet- 
ric monitoring will be presented in a forthcoming paper 
(Olsen et al. in prep) with other candidates followed. 



5. HD 168443: Superplanets in disks? 

The inner planet orbiting HD 168443 was detected by 
Marcy et al. (1999| ). They had gathered 30 HIRES/Keck 
radial-velocity measurements over 800 days that allowed 
them to characterize a 58-d period orbit from which they 
inferred a 5 Mj up companion to the star. These authors 
also mentioned a significant drift (89.4 ms- 1 yr~ 1 ) of the 
observed velocities indicating the presence of an addi- 
tional outer companion in the system, as yet undetected 
directly. Their careful check for a possible stellar compan- 
ion to HD 168443, in spectroscopic, adaptive optics and 
HIPPARCOS astrometric measurements, allowed them to 
constrain a potential stellar companion to be at a distance 
between 5 and 30 AU fro m the primary star and with a 



mass smaller than 0.5 M . Marcy et al. (1999| ) also men- 



tioned a clear indication of curvature in the radial- velocity 
drift that allowed them to postulate that the 2 nd compan- 
ions had to be on an orbit with a period of at least 4 years 
and be more massive than 15Mj up . 

After the announcement of the detection of a planet 
orbiting HD 168443, the star - part of our planet-search 
programme in the southern hemisphere - has been fol- 
lowed regularly with CORALIE. We gathered 58 additional 
observations over 670 days reaching the 2 nd extremum of 
the curve of radial- velocity variation due to the 2 nd com- 
panion. Combining our o wn measurements with the 30 
published Keck velocities (Marcy et al. 1999) we obtained 
the complete description of the system for the IAU 202 
symposium on Planetary Systems in the Universe where 
we presented the simultaneously-derived Keplerian solu- 



tions for the 2 companions (Udry et al. 20001). The sec- 
ond companion was found to be on a 1667-day orbit and 
have a m2 sini of 15.1 Mj up . This early solution is recalled 
in Table |. 

The star was of course closely followed by Marcy and 
collaborators as well. Within a few weeks they derived 
|a complete solution with their own data~( Marcy et al 
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Top. Simultaneous two-Keplerian solution for 
1684 43 system derive d from older Keck veloci- 



om Marcy et al. 1999 , open symbols) and sub- 
CORALIE measurements (filled circles). Bottom. 
around the combined solution 



duals 



2001) confirming our early result. On our side we also went 
on gathering 14 additional CORALIE observations in 195 
days. In Fig. we present an updated two-Keplerian com- 
plete solution of the system derived from the 72 CORALIE 



velocities and the early 30 Marcy et al. (1999) measure- 
ments. The derived parameters are given in Table |^. The 
quality of the fit is fair (x? e d = 2.46). The increase of the 
monito ring span since the first proposed solution allowed 
us and Marcy et al. (200T| ) to better cover the long-period 
orbit, estimat ed now to have a period of P = 1739 days 
(1770 days for [Marcy et al. 200l| ), slightly longer than the 
first estimate. Setting the primary mass to 1 M , the min- 
imum masses of the inner and outer companions inferred 
from our model are now 7.7 and 16.9 Mj up , respectively. 

5.1. Nature of the companion to HD 168443? 

The properties of HD 168443, if not due to unlikely orbital- 
plane inclinations for both companions, set very interest- 
ing qu estions on the nature and possible formation of such 



system s. Arc the companions of HD 168443 superylanets 
formed in the protoplanetary disk, brown dwarfs or even 
low-mass stars? 

The hierarchical organization of the system does not 
allow stability criteria to set constraining upper limits for 



the companions masses. Marcy et al. (2001) have shown 
that even a triple stellar system can be stable in the case of 
a coplanar geometry. Different orientations of the two or- 
bital planes favour, however, substellar companions. From 



Fig. 6. 

ments 
noise e 
best 
a linca 
Botton: 



Top. CORALIE temporal radial- velocity measure- 
}f HD 202206 (except observations with photon- 
■ror larger than 15 ms -1 ), superimposed on the 
rrjodel including a Keplerian orbital solution + 
radial-velocity drift with 42.9ms~ 1 yr _1 slope. 
Residuals around the combined solution 



(2001) also estimate the maximum mass of the outer com- 
ponent to be smaller than ~42Mj up . 

If formed outside the disk, the more massive outer com- 
panion is close enough to have truncated the disk within 
the ice limit, preventing thus a giant planet from forming 
in the outer regions and then moving towards the system 
center as predicted by the migration scenario. The inner 
companion has then also to be a brown dwarf. On the 
other hand, if the two companions are in the disk, they are 
located well within the ice limit of the young protoplan- 
etary disk. According to the migration scenario both ob- 
jects had thus to move simultaneously towards the central 
region of the system. Simulations by Kley (2000 ) explore 
this possibility and show that in such a case, the outer 
object accreates disk material more efficiently This rises 
the possibility of creating superplanets in the disk with 
masses larger than 15 Mj up . In such a case, if not very rare, 
these objects with masses above the high- mass tail of the 
observed planet-mass distribution (Fig. [l| |Jorissen et al 



2001) would represent a new population in the diagram. 



the HIPPARCOS astrometric measurements Marcy et al 



This seems to be rarely observed. The "final" answer will 
be given by future precise astrometric measurements that 
will determine the true masses of the companions. 

6. HD 202206: Triple system or superplanet in a 
binary? 

The CORALIE observations of HD 202206 started in 
August 1999. The obvious variation of the radial veloc- 
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Table 3. Same as Table g but for HD 168443. The early IAU 202 complete solution is given as well as the new updated 
one 



Parameter 



HD 168443 b 



HD 168443 c 



IAU 202 (Udry et al. 2000b) 



HD 168443 b HD 168443 c 

Updated orbit 



P 

T 
e 
V 
id 

K [ms" 

■At; mcs 

a(0-C) [ms- 

Xrcd 



[days] 

[JD-2400000] 
fkms" 1 ! 



58.117 ±0.014 1667 ±48 
51558.36 ±0.12 50269.5 ±36.0 
0.526 ±0.008 0.265 ±0.049 

-48.744 ± 0.002 
172.2 ±1.1 59 ±5 

473 ±6 288 ±13 

30" ± 58 b 
7.6 



58.116 ±0.001 1739.50 ±3.98 
51616.36 ±0.02 52014.5 ±3.6 
0.529 ±0.002 0.228 ±0.005 

-48.647 ±0.002 
171.61 ±0.22 63.67 ±0.84 
475.7 ±1.3 293.8 ±2.3 

30" ± 72 6 
3.0 a /8.1 b 
2.46 



ai sin i 
f(m) 
m2 sin i 
a 



[AU] 

[10~ 6 M e 

[M Jup ] 

[AU] 



0.002138 
0.386 
7.2 
0.29 



0.04559 
4.164 
15.1 
2.67 



0.002157 
0.396 
7.7 
0.29 



0.04574 
4.218 
16.9 
2.85 



1 HIRES data, 



CORALIE measurements 



ities allowed us to announce the detection of a low-mass 
companion to the star already after one orbital period, 
in May 2000 flESO 2000D , at the same time as HD 162020 
also described in this paper. When a second maximum of 



the ra cial-velocity curve was reached we noticed a slight 
drift of its value. We have now gathered 95 measurements 
covering more than 3 orbital periods. A simultaneous fit 
of a Keplerian model and a linear drift yields a period 
of 256 days and a large semi-amplitude of 565 ms -1 for 
the orbital solution (Table |^; Fig. 6). The orbit is fairly 
eccentric with e = 0.429. Choosing the stellar mass to be 
1.15 Mq (see above), the inferred minimum mass for the 
secondary is 17.5 Mj up . 

The quality of the solution is good with a weighted 
r.m.s. of 9.5 ms _1 around the fitted model (to be com- 
pared to the typical photon- noise error of 8 ms' 1 ) and a 
reduced x 2 value of 2.38. 

The slope of the radial-velocity drift is found to be 
42.88 ms'V" 1 - The available two older CORAVEL mea- 
surements obtained in 1989 and 1991 unfortunately do not 
allow us to further constrain the longer-period companion. 

Contrary to HD 110833 w hich was detected w ith a com- 
parable m 2 sini companion (Mayor et al. 1997 ) and then 



eccentricities involves the gravitational perturbation of a 
stellar companion, that can be applied to HD 202206. 

A more interesting characteristic of the system is given 
by the star's very high metallicity ([Fe/H] =0.37). (Bantoi 



et al. (2001a) have shown that a large fraction of stars with 



shown to be in reality a stellar binary ( |Halbwachs et al 



2000b|, the distance of HD 202206 (46.3 pc) prevents the 
HIPPARCOS astrometric data from constraining the vi- 
sual orbit. At such a distance the expected displacement 
on the sky of the star due to the inner companion is only 
0.26 mas, supposing sini = 1. A factor of 5 on the sin? 
bringing the companion into the stellar domain would still 
be insufficient at the HIPPARCOS precision. 

If not due to unfavourable orbital inclination, the ob- 
served low secondary mass sets the companion close to 
the limit of the planetary and brown-dwarf domains. The 
fairly large measured eccentricity does not help us to fur- 
ther constrain the nature of the object as most of the ex- 
trasolar "massive" planetary candidates are found on elon- 
gated orbits. An often proposed explanation for the planet 



this level of metal content have giant-planet companions. 
Once again, precise astrometric measurements will clarify 
the problem. 

7. Discussion 

Several scenarios may be invoked for the formation 
of companions to solar-type stars with masses in the 
planet/brown-dwarf transition domain. On the one hand, 
the observations of very light free-floating objects in young 
formation re gions ( Zapatcro Osorio et al. 200C ; Lucas fe 
Roche 200C) suggest that the minimum mass of brown 



dwarfs formed as stars by fragmentation of a protostellar 
cloud could be very small, at the ~ 5 Mj up level. On the 
other hand, several processes have been proposed for the 
formation of massive planets in stellar accretion disks: by 
gravitational instability of the disk (possibly triggered by 
the perturbation of an additional stellar companion, Boss 
[2000 , 2001) or by gas accretion of simultaneously migrat- 
ing planets trapped into resonances (e.g. Kley 2000). The 
interesting point in the latter scenarios is that a supple- 
mentary companion - stellar or planetary - is required or 
at least is supposed to enhance the process. 

7.1. The "companion's" influence 

To check if the potential companion's influence can be 
seen in the data, we made an inventory of the known 
visual companions of the planet-hosting stars. This vi- 
sual companion census was done through the Simbad-Vizir 
database (star within 1 arc-minute) and through the re- 
cent literature on adaptive-optics programmes searching 



for faint companions to stars bearing planets (Lloyd et al 
[2000 ; Luhman & Jayawardhana 2002). We also looked 
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tary mass distributions for stars hosting planets a. with- 
out further visual or spectroscopic companion (solid line), 
b. with a visual companion (dotted line), c. with an ad- 
ditional radial-velocity drift (long-dash line) unveiling a 
second planetary or stellar companion and d. with a vi- 
sual or a spectroscopic companion (i.e. b or c, short-dash 
line) . From these curves we directly see that the " spectro- 
scopic" and "visual" characteristics have the same effect. 
On the other hand, the curve for the stars without fur- 
ther companions rises less rapidly than the others, show- 
ing that "multiple" systems tend to harbour lighter plan- 
ets. This result is however not statistically significant. The 
Kolmogorov-Smirnov probability that curves a and d come 
from the same underlying population is 0.35. The curves 
even become indistinguishable if we restrict the sample to 
periods smaller than 1 year, so avoiding the bias favouring 
massive planets with long periods for which an additional 
radial-velocity drift is harder to see than for short-period 
systems, over the typical span of the current planet-search 
programmes. 



Fig. 7. Upper panel. Cumulative functions of the plane- 
tary mass distributions for stars hosting planets: a. with- 
out further visual or spectroscopic companion (solid line) , 
b. with a visual companion (within 1 arc-minute, dotted 
line) , c. with an additional radial- velocity drift (long-dash 
line) and d. with a visual or a spectroscopic companion (b 
or c, short-dash line) . Lower panel. Ratio of the number 
of stars hosting planets with visual or spectroscopic addi- 
tional companions to the number of stars without further 
companion, in 2-Mj up smoothing windows 



for indication of additional radial-velocity drifts in the 
known planetary orbital solutions, mentioned in the dis- 
covery papers or appearing in our CORALIE data. The 
census is probably far from being exhaustive but never- 
theless can bring initial insight into the question. The 
ratio of the number of stars hosting planets with visual 
or spectroscopic additional companions to the number of 
stars without further companions, in 2-Mj up smoothing 
windows, is shown in the lower panel of Fig. |^[ Unlike 
what was expected from the above-mentioned scenarios, 
there is no evidence that additional companions would 
favour massive-planet formation. There is even a trend for 
lower- mass planets (rr^sinz < 4Mj up , solid line) to ap- 
pear more often in "multiple" systems than massive ones 
(m-2 sini > 4Mj up , dashed line), by a factor of roughly 2. 
The statistics become too poor above 9Mj up to be able 
to draw conclusions about the more massive candidates. 

Is this observed difference significant? The upper panel 
of Fig. ^ presents the cumulative functions of the plane- 



7.2. Properties of mass subclasses 

It is worth noticing here that, in the lower panel of Fig. ^, 
the change in the mass regime is around 4Mj up , at the 
same position as the limit between the two possible modes 
of the planetary mass distribution in Fig. ^ (lower panel). 
This suggests that the bimodality of the mass distribu- 
tion is probably also not statistically significant and is 
potentially due to the mentioned observational bias. It is 
however interesting to check if the stars in the two mass 
subclasses show peculiar characteristics. This is achieved 
by comparing the distributions and cumulative functions 
of stellar and orbital properties for the two populations 
(limit at 4M Jup , Fig. |). 



7.2.1. Metallicity 

The metallicity distributions for the two mass subclasses 
are presented in the left column of Fig. |[ On average, the 
heavier companions seem to orbit stars slightly less metal 
rich than lighter planets (middle panel). The difference 
is however not statistically significant (KS prob = 0.061, 
lower panel) as already pointed out by [Santos et al. 
1(20014 ^). The same result holds when restricting the sam- 
ple to the CORALIE programme (KS prob = 0.065) and 
the same trend is also observed for the Lick survey by 



9 In order to decrease the statistical noise, the distributions 
in Figs. ^ and |^ are shown for all detected exoplanets. It has 
however been verified that the results do not change qualita- 
tively when considering only the stars in the volume-limited 
CORALIE planet-search sample (41/78 stars) 



Fischer et al. (2002) who compare the mean metallicities 
of two mass classes (limit at 5Mj up ). The difference origi- 
nates mainly from a decrease with increasing metallicities 
of the distribution of stars with "massive" companions 
whereas the distribution for stars with lighter planets in- 
creases in the same metallicity range. An improvement of 
the available statistics will shed light on the question. 

An interesting feature of the metallicity-mass distri- 
bution is the lack of massive planets at very low metal- 
licities (upper left panel), although stars with massive 
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Fig. 8. Metallicity, eccentricity and period distributions of star hosting planets. Comparison between subsamples with 
different planetary masses: mi sini < 4Mj up (filled circles, open histogram and solid line) and m-i smi > 4Mj up (open 
circles, filled histogram and dotted line). The Kolmogorov-Smirnov probability that the two distributions come from 
the same underlying population is given in the lower panels. The probable brown dwarfs HD 114762 (in the upper- left 
panel) and HD 162020 (in the upper-right panel) are represented by crosses superimposed to open circles 



planets are on average more metal poor than stars with 
light planets, as seen above. Below [Fe/H] ~ —0.25, only 
HD 114762 (7712 sini = 11 Mj up ), often considered a brown 



7.2.2. Eccentricity-Period 



dwarf ( j(Jochran 1991 j, has a mass larger than 1.1 Mj up {$ 
in the upper lett panel ot fig. ^j. Un tire low-metaliicity 
side of the diagram, there even seems to be a limit im- 
posing a minimum metallicity for the star to harbour a 
planet with a minimum given mass. This corroborates the 
idea that more "solid" material is needed in the accretion 
disk to form more massive planets. It also may be inter- 
preted in terms of shorter time scales needed to accrete the 
planet core in metal rich environments, leaving more time 
for the planet growth over the lifetime of the disk. Both 
interpretations support the gas-accretion scenario for the 
formation of giant planets. 



It has often been pointed out (e.g. Heacox 199E ; Mayor & 



Udry 2000; Udry et al. 2000b; Btcpinski fc Black 2001 , and 



references therein) that, for periods larger than a few tens 
of days, the eccentricity distributions of planetary systems 
and stellar binaries are unexpectedly similar. For shorter 
periods, evolutionary effects (planet migration, tidal cir- 
cularization) change the distribution, favouring low- or 
zero-eccentricity orbits for both populations. The middle 
column of Fig. || shows that this is also true for the two 
defined subclasses of planetary masses. The distributions 
presented in the middle panel are very similar except for a 
prominent peak of "light" planets at small eccentricities. 
The latter corresponds to short-period, close-in planets, 
probably circularized through the migration process. They 
also relate to the accumulation at short periods in the mid- 
dle panel of the 3 rd column in the figure. When restricting 
the sample to periods longer than 50 days, the mentioned 
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peaks - mainly responsible for the difference observed in 
the cumulative functions (KS-prob values in the figures)- 
disappear. Note however that the longer-period part of the 
distribution is still very much observationally biased. 

An important point to note here is that massive plan- 
ets on short-period orbits are rare - although they are 
more easily detected than lighter ones. If we consider the 
companion of HD 162020 ((g> in the upper right panel) to 
be a brown dwarf as shown in Sect. 4, the only observed, 
short-period massive candidate is r Boo, close to the cho- 
sen limit of the considered mass subclasses. While sub- 



mittin g this pap er, we learned about a study by |Zucker 
& Maz jch (20021) estimating the statistical significance of 
this feature. 



7.2.3. Concluding remark 

We also have searched for differences between other dis- 
tributions of orbital and stellar-host properties of "light" 
and "massive" exoplanets, without success. In particular, 
the primary masses do not correlate at all with the planet 
masses. The only observed trends are the need for metal- 
rich stars and the lack of short periods for massive planets. 
If confirmed with improved statistics, these features may 
bring constraints for the migration scenario. Possible ex- 
planations may invoke the idea that massive planets do 
not migrate as easily as lighter ones or, on the contrary, 
that they cannot stop their migration process when reach- 
ing the central part of the system, falling into the star. 
The higher metallicity of stars hosting light planets may 
support this latter view. 
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